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QUANTUM STATISTICAL CORRECTIONS TO ASTROPHYSICAL PHOTODISINTEGRATION RATES 
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ABSTRACT 

Tabulated rates for astrophysical photodisintegration reactions make use of Boltzmann statistics for the pho- 
tons involved as well as the interacting nuclei. Here we derive analytic corrections for the Planck-spectrum 
quantum statistics of the photon energy distribution. These corrections can be deduced directly from the 
detailed-balance condition without the assumption of equilibrium as long as the photons are represented by 
a Planck spectrum. Moreover we show that these corrections affect not only the photodisintegration rates but 
also modify the conditions of nuclear statistical equilibrium as represented in the Saha equation. We deduce 
new analytic corrections to the classical Maxwell-Boltzmann statistics which can easily be added to the re- 
verse reaction rates of existing reaction network tabulations. We show that the effects of quantum statistics, 
though generally quite small, always tend to speed up photodisintegration rates and are largest for nuclei and 
environments for which Q/kT ^1. As an illustration, we examine possible effects of these corrections on the 
r-process, the rp-process, explosive silicon burning, the 7-process and big bang nucleosynthesis. We find that 
in most cases one is quite justified in neglecting these corrections. The correction is largest for reactions near 
the drip line for an r-process with very high neutron density, or an r/j-process at high-temperature. 

Subject headings: Abundances, nuclear reactions, nucleosynthesis 



1. INTRODUCTION 

The capture of a projectile particle by a nucleus followed 
by the emission of a photon is called radiative capture. The 
inverse process is called photodisintegration. Both types of re- 
actions play important ro l es in stellar and big ban g nucleosyn- 
thesis dRolfs & Barnes 1990; Smith, Kaw 
19931: Waaoner, Fowler, & Hoyl^ 

Fowler, Caughlan, & Zimmerman 19671 
Iliadisi.2007,) . Photodisintegration in astrophysical environ 





1967; 


IClavtonI 


1968; 



ments often involves a thermal distribution of photons that 
excite nuclei above the particle emission threshold. Above a 
temperature of T ^ lO** K, photodissociation can become a 
dominant process in the reaction flow and the photo-ejected 
nucleons can be captured by other nuclei leading to photodis- 
integration rearrangement as can happen for example during 
core or explosive oxygen or silicon burning (IClaytonlll968t 
ETiadis 2007). 

In determining the rate for photodisintegration reactions, 
however, one should take into account various factors aris- 
ing from dealing with photons in a two body problem. One 
of them is the fact that photons are massless bosons, and 
hence, obey Planckian statistics. Historically, photodisinte- 
gration reactions have been treated with Maxwell-Boltzmann 
statistics, bo th because it is usually an excellent approxi - 
mation (cf. iRauscher, Thielemann & Oberhummeil (119951) ; 
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IlliadisI ( 12007 1)) and because this assumption simplifies the de- 
termination of the photodisintegration rates. However, low- 
energy photons are more correctly represented by a Planck 
distribution. As discussed below, if the ratio of the capture 
2-value to the temperature Q/kT is large then the use of a 
Maxwell-Boltzmann distribution for the photons is a good ap- 
proximation. On the other hand, for Q/kT ^ 1 (e.g., if one 
is interested in the photo-ejection of loosely bound particles 
from nuclei), then the effects of quantum statistics become 
more relevant. We show here that modified thermal pho- 
todisintegration rates can be written in the form of a small 
analytic correction to the tabulated reaction rates obtained 
with Maxwell-Boltzmann statistics. Furthermore, these cor- 
rection factors are nearly independent of the nuclear cross 
sections and to leading order only depend upon the reaction 
Q-value, the Gamow energy (for the charged-particle nonres- 
onant part), and the resonance energy (for the resonant part). 

1.1. Photonuclear Reactions 
For a reaction (not necessarily in equilibrium) of the form 



1+2^3 + 7 , 3 + 7^ 1+2 , 

the forward and reverse reaction rates are given by 

ri2 = nin2{<Jv) 12 = niXn 
rj2 = n3n^{ac)ji = nT,Xj2 . 



(1) 



(2) 
(3) 



Here, (crv) 12 denotes the thermally averaged reaction rate per 
particle pair for the capture reaction and is given by an inte- 
gration over an appropriate velocity distribution 0(v), 



crv) 12 = / 0-12 (v) v(f>(v)d V , 



(4) 



where v denotes the relative velocity of the nuclei 1 and 2. For 
most cases of interest in astrophysics, the massive interacting 
particles are non-degenerate (i.e., dilute) and non-relativistic 
(i.e., their rest mass energy is large compared to kT). Hence, 
one can use a Maxwell-Boltzmann velocity distribution with 



2 



Newtonian kinetic energy because 
1 



\J P^c-+m^c^ jkT ^g-'jk 



(5) 



where the exponential factor involving the rest mass energy 
drops out with a proper normalization of the distribution 
function. Also note that when both particles 1 and 2 obey 
Maxw ell-Boltzmann statistics, so does their relative velocity 
(IClayton.1968) . Hence, the thermally averaged capture rate 
per particle pair is given by 



{av)n=\l — (kT)-i 

TT/i 



ai2{E)e-'^l" EdE 



(6) 



where ji denotes the reduced mass for the particles 1 and 2. 

In the photodisintegration reaction rate, however, the rela- 
tive velocity of the photon with respect to the target nucleus is 
always the speed of light c. This eliminates any dependence 
of the react ion rate on th e velocity distribution of the target 
nuclei (Th ielemann et alj |1998). Therefore, the photonuclear 
reaction rate {ac).yj, becomes an integral of the reaction cross 
section over a Planck energy distribution for the photons: 



1 



2Ci3)ikTy Jmax(0,G) 



CT^3 (£7) 



1 



■ Et dE^ 



(7) 

Using the fact that for a Planck distribution 



fkT 

we can equivalently write Eq. (|7]i in more familiar form 

stt n ^ 1 



(8) 



A 



73 ■■ 



(hey 



max(0,e) 



(T73 (£7 



E]dE-^ , (9) 



where C(3) = 1.20206 is the Riemann zeta function and E-y 
denotes the photon energy. The integration threshold is the 
2-value of the capture reaction (see Figure [U or zero in the 
case of negative Q. 

1.2. Detailed Balance Condition 

It is difficult to determine the cross section 0-73 directly 
from experiment. However, the interaction between photons 
and matter is very weak {e^ /Tic ^ 1) so that the reaction can 
be treated with first order perturbation theory. In this case, 
the transition probabilities become proportional to the matrix 
elements of the perturbing Hamiltonian and the hermiticity of 
the perturbing Hamiltonian gives rise to a simple relation be- 
tween the capture and disintegration cross secti ons. This is 
known as the detailed balance equation (B latt & Weisskop^ 
1199 ll) . For a reaction involving a ground-state to ground-state 
transition for two nuclei with energy E, leading to a gamma 
ray with energy E^=E + Q this is given by 



(7^i(E^) = 



gigi ^ic^E 
g^(\ + 5n) E^ 



cJiiiE) 



(10) 



where g, = 2 j, + 1 are the spin degeneracy factors for the 
ground state of the nuclei and the Kronecker delta function 
accounts for the special case of indistinguishable interacting 
nuclei. Using this detailed balance equation, the photodisin- 
tegration rate for a single-state transition can be related to the 
forward capture rate. 



Substituting Eq. ( fTOt into Eq. (|7]) and also changing the 
variable from £7 to £ = £7 - 2 in the integration, one can 
average over the velocity distribution of the ground-state in- 
teracting nuclei {ac)^3 as follows: 



{0-C)j3 ■■ 



8182 



<yi2(E)- 



1 



-EdE. 



(11) 



2a3)(kTyg,(l+Sn)Jo ^e(£+G)Ar-i 
At this point, one usually introduces the approximation: 



Here, we point out that by inserting this approximation and 
then correcting for it, Eq. ( fTTI ) can be rewritten in the follow- 
ing exact form: 



(CTC 



\ \ f v2^\ glgl ( 



3/2 



-Q/kT 



(13) 



where 7? is a small and dimensionless number which is for- 
mally given by 



\+R: 



J^^an(E)(e'-^^Q^/''^-ir'EdE 



J^^ai2{E)e-(E+Q)/kTEdE 



(14) 



1.3. Thermal Population of Excited States 

The generalization of Eq. ( fTJt to the average over ther- 
mally populated states among the initial an d final nuclei is 
straightforward dClavtonll 19681: Iniadisll2007[) . One must first 
replace the ground state (g.s.) to g.s. forward reaction cross 
section <ti2 with a weighted average over the thermal popula- 
tion of states /i in the target nucleus 1, and also sum over all 
final states in product nucleus 3. (Note that we only consider 
light particle (p,j), (n,7), or (0,7), reactions for which we 
can ignore their excitation.) Thus, the effective stellar ther- 
mal forward rate becomes 



-Ei^JkT 



E 



(15) 



which can also be written as 



Mt2 = /?«Mi2 , (16) 

where the stellar enhancement factor R„ is defined by 



..-EuJkT E.(o-i2i')r2 



-E,^/kT 



(17) 



Usually, tabulated thermonuclear reaction rates are given 
as the ground state rate and the stellar enhancement fac- 
tor must be determined from a statistical model calcu- 
lation as in Holmes et al. (1976), Wooslev et.al. (1978), 
iRaus cher & Thielemann (.2000) and .Rauscher & ThielemannI 
(2004. 

The thermally averaged photonuclear rate for a distribution 
of excited states in initial and final heavy nuclei then becomes 



crc 3.V = ■ 



Y.^831' 



-Ei„lkT 



I/— >-/J, 



(18) 



The generalization of the detailed balance condition of 



3 



Eq. O: 



kT 



(19) 



where 7?^^ denotes the use of Q^i, and a'^^^{E) in Eq. ( fT4l l. 
Inserting Eq. ( fT9] l into Eq. (fTSl l and using the fact that = 
Q-Ei^ + Eif^, we can write 



(ac);3=(l+/?)(av)t2 



8C(3) 



3/2 



G\G2 



G3(l+<5l2) 



-e/*7- 



(20) 



where /? represents an average correction factor among all 
thermally populated states. As demonstrated below, R is 
nearly independent of the detailed nuclear structure. Hence, 
we can simply utilize the ground-state Q-value as a represen- 
tative average over the distribution of Q-values among the 
thermally populated states. Also, the spin factors above are 
now replaced by the relevant nuclear partition functions G, : 



IkT 



(21) 



where a denotes the individual states in nucleus /. 

Stellar reaction rate tables are usually listed as func- 
tions of temperature Tg in units of 10'' K and are given as 
[NA{(J\iT>)))*] . Thus, we can rewrite Eq. ^ as 



(22) 



X9.8685 X 10'' (/trg) 



3/2. 



G^{l + 5n) 



-ii.eosQ/r, 



where now A'a is Avagadro's number so that [A'A(crv(79))*] is 
in units of cm-' mol~' s~', Q is in units of MeV and jl is the 
reduced mass in atomic mass units. 

Equations ( |20l ) and (|22] ) are in a convenient form 
because in the limit of 7? — > 0, they reduce to usual 
photodisintegration rates available from various com- 
pilati ons (e.g.. Fowler, CaughlarK_& Zimmerman (1967, 
1975h: [H olmes et al. (197ii^ |m)oslev et.al. (1978); 
Caughlan & Fowler (1988); NACRE7Angulo et al. (1999); 
TALYS, iGo rielv. Hilaire, & Konig (2008); NONSMOKER, 
iRauscher & Thielemanra (,2000) or REACLIB, Cvbu rt et alj 
( 120101) ). The combined factors multiplying (crv) *2 in Eq. ( [22] l 
are usually referred to as the "reverse ratio" as this factor 
gives the reverse reaction rate in terms of the forward rate. 
In this work, we show that there is a simple correction 
{l+R) to this reverse ratio due to the difference between 
Planckian and Maxwell-Boltzmann statistics. For most of 
the remainder of this manuscript, our goal will be to derive 
a simple analytic form for R for ease in correcting existing 
tabularized reverse reaction rates. We will also derive simple 
analytic approximations to clarify the essential physics of this 
correction and summarize examples of which astrophysical 
conditions may be most affected by these correction factors. 

1.4. Nuclear Statistical Equilibrium 

Before leaving this discussion, however, it is worth empha- 
sizing again that the above rate does not imply equilibrium. 



but only detailed balance and a thermal population of pho- 
tons and nuclear excited states. Nevertheless, the situation of 
equilibrium between capture and photodissociation frequently 
occurs in astrophysical environments and is referred to as nu- 
clear statistical equilibrium (NSE). It is of note that the condi- 
tions of NSE are also modified from the usual Saha equation 
by the above quantum corrections. Moreover, in conditions of 
NSE, one sometimes synthesizes nuclei for which Q/kT ^ 1 
and the corrections can become larger. Examples of this in- 
clude the formation of nuclei near the proton drip line in the 
hot hydrogen burning rp-process, or the synthesis of nuclei 
near the neutron drip line in the neutron-capture r-process, as 
discussed below. 

To see the revised conditions of NSE consider the evolution 
of a nucleus undergoing rapid particle captures and photodis- 
sociation. This can be written as 



dt 



-n 1 n2 (ctv) 12 + n3n7 (o-c) ^3 



(23) 



The equilibrium condition {dn i/dt) = therefore demands 
that 



nin2 n^K'^cifi 



«3 



/I _L^.^/^27^/^^^y''^ G1G2 



-Q/kT 



(24) 



or in terms of mass fractions and temperature, it is more con- 
venient for stellar models to write 



±}±l = 9M85xlO'(l+R)^^^ 
X3 ^ ^ p 



G1G2 



G3(l+(5l2) 



-I1.605Q/r9 



(25) 



In the limit that R ^ 0, Eqs. (|24li and represent the 
usual nuclear Saha equation ( Saha, 192 1) of statistical equi- 
librium which also invokes the Maxwellian approximation 
given in Eq. ( fT2] ) either directly or indirectly in its derivation 
(cf. Clayton ( 1968); Iliadis (2007)). The deviation of NSE due 
to quantum statistics may impact the evolution of explosive 
nucleosynthesis environments for which one can encounter 
nuclei with small photodissociation thresholds, e.g., near the 
neutron or proton drip lines. To the extent that such nuclei are 
beta-decay waiting points, for example, the altered statistics 
will affect the timescale for the build up of abundances. An- 
other possible application of the corrections deduced here is 
for the ionization equilibrium of atomic or molecular species 
with a low ionization potential in stellar atmospheres. How- 
ever, we will not consider that case further here. 

2. EVALUATION OF THE CORRECTION FACTOR R 

It is worth noting that quantum effects always tend to speed 
up photodisintegration rates because the Planck distribution 
places many more photons at low energy than a Maxwell- 
Boltzmann distribution of the same temperature or energy 
density. In other words, R is always positive definite. Hence, 
even though it is often small, it is worth including. This cor- 
rection could of course always be evaluated by direct numeri- 
cal integration of Eq. ( fT4l i or Eq. ( fTTI) . In a large network cal- 
culation with evolving temperature, however, the repeated nu- 
merical integrations would slow the computation time, more- 
over, it is a tedious task to assemble all of the relevant cross 
section data. 

Nevertheless, the advantage of introducing the R in Eq. ( |20] | 
is that it gives the quantum correction as a small fraction of 
the classical Maxwellian result. Hence, for implementation in 



4 



large networks, an analytic approximation to the exact numer- 
ical integration for R is both adequate and desirable. More- 
over, we show that an accurate analytic correction is readily 
available based upon the input from existing reaction rate ta- 
bles (either in analytic or tabularized form). We also show that 
R is nearly independent of the nuclear cross sections, and to 
leading order only depends upon the value and the Gamow 
and/or resonance energy. 

The key to evaluating R is to perform a geometric series 
expansion of the numerator of Eq. (fl4l l. i.e.. 



^ oo 

J— = ye-"\ 



(26) 



temperature decreases, as well as the n"-*/^ pre-factor and the 
exp[-(n- 1)2/^7^]- These three conditions guarantee that 
this is a well behaved, rapidly decreasing, convergent series 
even as g ^ 0. In what follows, we show some illustrations of 
the magnitude of these corrections and also derive some alter- 
native analytic forms to illustrate the basic structure of these 
corrections in more detail. Indeed, we show that in practice, 
only one or two terms are needed in the series and the correc- 
tion factors are largely independent of the underlying nuclear 
structure. As a caveat to the reader, however, we note that 
the analytic approximations derived below do not include the 
stellar enhancement factors, and as such should be used with 
caution in a real astrophysical plasma. 



Note that terminating this series after the first term corre- 
sponds to the usual Maxwell-Boltzmann approximation in 
which case one obtains R = as expected. Using the full ex- 
pansion, however, leads to 



11=2 



(27) 



(28) 



where the coefficients r„ (for n > 2) are given by 

(j*{E)e-"'^l"EdE 

Eq. dZTl i gives /? as a geometric series in powers of e~^l''^ 
which is usually a small quantity. The coefficients r„ are also 
less than unity and they rapidly decrease with increasing n. 
Hence, the series given in Eq. ( |27] | is a rapidly convergent 
one. However, the explicit determination of the r„ requires 
some attention to the energy dependence of the cross sec- 
tions. Nevertheless, this task is greatly simplified when the 
forward thermonuclear reaction rates as a function of temper- 
ature {(jv(kT)) have already been compiled. Combining the 
expression for the reaction rate in Eq. (|6]l with Eq. (l28T l. the 
expansion coefficients become 



1 {av{kT/n)Y 
''"'1^ {av{kT)Y 



(29) 



This immediately gives the correction factors in terms of the 
tabulated rates. 



R- 



E 

n=2 



]_ (MI£My-(„-l)Q/kT 



(3/2 {av(kT)Y 



(30) 



or converting to units of and noting that reaction rate com- 
pilations are in terms of [A^A(fv(79))*] we have 



R- 



1 [NA{(Jv{T^ln)Y\ 
^n3/2 [NA{av{Tq)Y] 

n=2 



-I1.605(n-I)e/ri, 



(31) 



These expressions for the reverse rate correction factor R 
make clear the physics of the correction factor The factors 
become a sum of the correction factors in terms of thermonu- 
clear averages over decreasing temperatures, kT — > kT/n. 
These terms achieve the task of increasing the photodissoci- 
ation rate due to the fact that the Planck distribution includes 
many more low energy photons as (E + Q)/kT 0, which 
is what the [e\p{{E + Q)/kT}- 1] denominator in Eq. ( fTTT i 
achieves. 

Eqs. ( l30l l and dSTI ) are the key equations for this paper. 
Note that {av(kT)) is a rapidly decreasing function as the 



2.1. Non-resonant Charged-particle Capture Reactions 

When the projectile particle is charged, the capture reac- 
tion must tunnel through the Coulomb barrier at low energy. 
To account for this, t he cross section can be f actored into the 
following form (Fowl er. Caughlan. & Zimmerman 1967 ): 



cr(£) = 5(£)^e" 



(32) 



Here, S{E) is called the astrophysical 5-factor and Eg is the 
Gamow energy which characterizes the penetrabihty: 



1 /27rZiZ2e2 



: 0.97907 W (MeV) , 



where 



w = z\zIp 



(33) 



(34) 



and fl is in atomic mass units. 

The S factor contains information about the detailed nuclear 
interaction. Away from resonances, the astrophysical S(E) 
factor is a slowly varying function, and at low temperatures 
the integrand is dominated by a small region known as the 
Gamow window (iFowler. Caughlan. & ZimmermanI 19671) lo- 
cated at an energy £0 as defined below. As such, to the 
desired accuracy it can be replaced with an average effec- 
tive value, 5'eff(£'o)- Therefore, at low temperature S(E) will 
cancel in the ratio given in Eq. ( |28] ). However, for higher 
temperatures the variation of the S'-factor with energy over 
the Gamow window can become relevan t. Hence, following 
iFowler. Caughlan. & ZimmermanI (Il967h we write 



S{E) = S{Q) 



i.'Me.'JMe^,. 

5(0) 2 5(0) 



(35) 



(36) 



where the dot denotes derivative with respect to energy. 

Inserting this cross section into Eq. ( l28b the correction co- 
efficient r„ becomes 

where the function F is defined as follows: 

poo 

F (Ec/kT) = / 5(£)e"V^e-^/*^t/£. (37) 
Jo 

This is a familiar integral i n nuc lear astrophysics 
( IFowler. Caughlan. & ZimmermanI 119671) . Eq. ( IJTT i 

corresponds to a product of the barrier penetrabihty 



5 



times a Maxwellian distribution. It is strongly maxi- 
m ized in the Gamow window and w e ll appro xi mated 

fowler. C aughlan. & ZimmermanI 119671: |Clayton' 119681: 
adis 2007^ as a Gaussian integral near the maximum of the 
integrand. Hence, we write 



S'-factor we have 



F(Ec/kT) = 



1/2 



A5eff(^r) e 



-iEo/kT 



where 



Eo = (VE^kT/2f/' 



= 0.12204 w'/^t;^''^ 

is the peak of the Gamow window and 

A = 4^(EokT/3) 

= 0.23682 W'/'To^'''^ MeV 



(38) 



(39) 



(40) 



is the peak width. The eff ective 5-factor for charged-particle 
reactio ns is given by iFowler. Caughlan. & ZimmermanI 
([T967h : 



S^flikT) = S(0) 

1 m 



l + ^ + ^{Eo + (35/36)kT) 



+ - ^(£o' + (89/36) 



(41) 



where t = 3Ei)/{kT). Converting to temperature, this be- 
comes 



5eff(79) = 5(0) 



l+0.09807rJ''V"'''^ 



+ 0.08378—^9 + 7.447 x IQ-^^W^^^T^^^ 
S S ^ 

+ O.OBOO^W'/Vg/'' MeV-barn. 



+ 0.1220^W>/^7g^/' 
;^(0), 



(42) 



The evaluation of the remaining terms in the series can 
be done simply by making the replacement kT ^ kT/n in 
Eq. (in. 

Now collecting all terms, the coefficient r„ becomes 

^^^S fT/n) ^_[3(../4..)---(„v-M)] ^43) 
n / •Seff(fcr) 

jSeff(V^^- [4.2487(W/r,)'/^ 

which immediately leads to the desired analytic expression 
for/?: 

^n^f'S^sikT) 



n=2 



E SsnjTg/n) -[ii.605(n-l)e/r9+4.2487(M'/r9)'''^(«'''^-l)] 



n=2 



where Q is in MeV in the latter equation. Note that the de- 
pendence on 5(0) cancels in the r„ and the correction factor 
R. Also note that as the ratio of the S'gff factors approaches 
unity, the correction factor only depends upon the Gamow en- 
ergy and the reaction Q-value. Hence, for a slowly-varying 



,J_^-[3(£./4.r)"-X«''^-i)] 

1 -[4.2487(H'/r,)'/'(n'/^-l)] 



(45) 



and 



R = \^ —e~ [(«-l)GAr+3 {Ea/4kT) (n"'-l)] 
„5/6 



«=2 



(46) 



E l -[ll.605(«-l)e/J9+4.2487(H'/r9)'''-^(«'''-^-l)] 
„5/6^ 



„5/6 
«=2 

In practice, this is a rapidly converging series requiring at 
most the first few terms for sufficient accuracy. To illustrate 
the rate of convergence, the first few r„ coefficients of Eq. i45[ 
are plotted as a function of Ec/kT in Figure |2] It is clearly 
seen in this figure that these coefficients become small as n 
increases and that the terms in the series are nearly of the 
same functional form. Hence, if speed is desired, an ade- 
quate approximation can usually be obtained by retaining only 
the first term in the series for R. For illustration the dashed 
lines in Figure |2] show a comparison of the exact integration 
of Eq. ( [36] l with the analytic expansion (Eq. ( |45] l) in the case 
of a slowly varying 5-factor. Clearly, the analytic series is 
adequate for the r„. The slight difference between these two 
sets of curves relates to the fact that the Gamow window is 
not exactly Gaussian even in the case of a constant 5-factor 
This, for example, is the reason that for a constant S factor as 
Ec —7- 0, the exact expression in Eq. (|36] l goes to 1 /n as can 
be seen from Eq. (|37| i. whereas the analytic approximation 
expansion goes to l/n^^^. 

Figure |3] illustrates the charged-particle reverse-reaction 
correction factor R from Eq. ( |46] | as a function of Q/kT for 
several values of the dimensionless quantity Ec / Q- The solid 
line is based upon the numerical integration of Eq. (fT4l l with 
a constant S(E) factor while the dashed line is based upon 
the first three terms of the approximate series given in Eq. 
(|46] |. Even though the r„ remain large up to Ec/kT > 10, 
the total correction R is only ^1% for reactions for which 
Q I kT^ 1 . From this figure it is also clear that retaining only 
the first three terms is an excellent approximation down to 
Q/kT ^ 0.4. Below that, however one should probably in- 
clude more terms in the expansion depending upon the Ec/Q 
ratio. We have found, however, that one never needs more 
than about six terms in the series even in the limit that Q^Q. 

As illustrations of corrections to the reverse ratio for an 
evaluated table, the plots on Figure |4] show the R factors for 
the hydrogen-burning charged-particle reactions ^Be(/:i,7)**B 
{Q = 0.135 MeV) and ''He(a,7)^Be {Q = 1.587 MeV) as a 
function of Ty. These plots were generated from Eq. ( [3T| l 
using the REACLIB compilation (Cyburtetal. 2010) for 
[NA{o'v){kT)] and reaction Q-values. Although most hydro- 
gen burning takes place at relatively low temperature (Tg < 1), 
evaluated thermonuclear reaction rates are usually tabularized 
up to Ti)= 10. Hence, we consider the same range here. 

On each plot, the solid line shows the application of 
Eq. dSTT i. the dashed line shows the application of Eqs. (|44] |. 
while the dotted line shows the results obtained from keeping 
only the first term in the series in Eq. ( |46] |. Surprisingly, the 
correction for the reverse ^He(a,7)^Be reaction based upon 
the REACLIB evaluation is almost indistinguishable from the 
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approximation based upon only the first term in the series with 
a constant S'-factor. This is due to the fact that the S factor in 
the compilation is taken to be constant at high energies and 
that Q/kT > I for the plotted range of Tg. These figures show 
that the corrections are generally quite small and adequately 
represented by only the first term in Eq. ( |46] |. 

2.2. Non-resonant Neutron-capture Reactions 

At low energies where the de Broglie wavelength of the 
neutron is much larger than the radius of the target nucleus, 
the non-resonant neutron-capture cross section is proportional 
to 1 /v. Hence, 



(e)Ve' 



! constant . (47) 

At high energies, however, deviations from Eq. ( |47] | occur 
Such deviations can be described in terms of a Maclaurin se- 
ries in E'/^ (Fowler, Caughlan, & Zimmerman 1967) which 
roughly accounts for the contribution of higher partial waves 
to the cross section: 



(Tl2V=S(0) 



SiO) 



2 5(0) 



(48) 



where the dot denotes a derivative with respect to Zi'/^. The 
integration of Eq. ( l48T l over a Maxwellian velocity distribution 
gives (.Fowler. Caughlan. & Z immerman . 1967b 



((Ti2v) = SefiikT) 
where S^fdhT) is defined by 



(49) 



S,s(kT) = S(0) 



1 + 



2 5(0) 



(kT) 



1/2. 



3 5(0) 
45(0) 



kT-\-. 



(50) 



= 5(0) 



5(0) , /, 5(0) 
1 +0.3312^T^(r9)'/2 + 0.06463--^r9 + . 



5(0)' ' ' 5(0) 
Inserting this into Eq. ( |27l ) then immediately gives, 

1 S^fs(kT/n) 



R = 



jf^nV^ S,ff(kT) 

oo 

y^g-ii.6Q5(»-i)e/r9 

n=2 



-(n-\)QlkT 



cm"'s ' 



(51) 



1 [l+0.3312|g^(r9/«)'/V0.06463g|g(r9/«)] 

n^'^ [l+0.3312|®(r9)i/2 + o.06463|®r9] 

in which the leading 5(0) factor obviously cancels. In the 
astrophysical environments of most relevance for the present 
application one can often ignore the derivatives of 5(0). We 
are also mainly concerned with nuclei and environments with 
Q/kT^\ for which Eq. ( BTT i quickly converges. Hence, re- 
taining only the first few terms in the expansion we can write 
a sufficiently accurate correction as 



R'r 



e 



-Q/kT 



e 



-IQ/kT 



■ + - 



e 



-iQlkT 



■ + - 



23/2 33/2 43/2 

g-11.605Q/r, g-22.210G/r>, g-34.815e/r9 



23/2 



■ + - 



33/2 



43/2 



(52) 



Figure |5] shows the reverse rate correction factor for nonres- 
onant neutron capture as a function of Q/kT . The solid line 
is for an exact numerical integration of Eq. (O. The dashed 



line is from the analytic expression given in Eq. (1521 ) trun- 
cated after the first three terms. The first few terms in Eq. ( |52] | 
are an adequate approximation until Q/kT^Q.2, below which 
one should include more terms. 

2.3. Resonant Capture Reactions 

Often one encounters charged-particle and neutron-capture 
reactions in which the thermonuclear reaction rates can be 
dominated by single (or a few) low-lying resonances. In such 
cases, (T12 in Eq. ( fT4b is replaced by the Breit-Wigner resonant 
capture cross section for each resonance. 



(Ji2{E) = 7rX- 



(E-Er)^ + (Ti,/4) 



(53) 



where A is the de Broglie wavelength, = gri^ + Si2)/igig2) 
includes spin factors of the reaction, while F, is the particle 
(e.g., proton or neutron) width, F-y is the width for gamma de- 
cay from the resonant state, and is the observed resonance 
energy. 

It is worthwhile to consider the limit in which the total res- 
onance width is small and dominates the reaction rate. In that 
case, in the limit of Ftot 0, we have 



1 



rtot/2 



7r(£-£,)2 + (F2,j/4) 



S{E-Er) 



(54) 



Inserting this into Eq. 
plified and reduce to 



for r„, the integrals are greatly sim- 



-(n-l)E,./kT 



(55) 



After summing the series, this leads to the final correction 
factor for a single resonance of 



R-- 



g(Er+Q)lkT 
g(E,+Q)lkT _ I ■ 



1 



(56) 



e 



11.6O5(£,.+0/r9 



gll.605(£,+G)/r5_l 



1 



where, in the second equation the resonance energy and Q- 
value are in units of MeV. As in the above cases, this correc- 
tion factor is only > 1 % in reactions for which {Er + Q) /kT ^ 
1. 

As an illustration of a resonant reaction. Figure |6] shows 
the reverse rate correction factor for the resonant reaction 
'2C(/?,7)13N {Q = 1.944 MeV) as a function of T9. The solid 
line on this plot was generated from Eq. (|3TI ) using the REA- 
CLIB compilation (Cyburtetal. 2010) for lNA{(7v){Tg)] and 
the reaction Q-\alue, while the dashed line shows the result 
from an application of the simple single resonance correction 
in Eq. ( |56] |. Even though this reaction has a second resonance 
at higher energy, most of the correction factor is accounted for 
by the single resonance approximation. 

3. APPLICATIONS 

Having deduced the analytic corrections it is worthwhile to 
briefly consider some illustrations of the practical applications 
of the above corrections. From the discussion above it is clear 
that these correction factors arise as Q/kT ^ I. For practi- 
cal applications in astrophysics, this implies T() -^l. We now 
consider several astrophysical environments in which these 
quantum corrections might appear These include the rapid 
neutron-capture reactions near the neutron drip line during 
the r-process, the rapid proton capture near the proton drip 



7 



line during the r/?-process, core or explosive oxygen or silicon 
burning in massive stars, the 7-process formation of proton- 
rich nuclei, and the first few moments of cosmic expansion 
during the epoch of big bang nucleosynthesis. 

3.1. r-Process 

The r-process involves a sequence of rapid neutron cap- 
tures in an explosive environment (" Burbidge et al.l Il957t 
lMathews~&" Ward 1985). It is responsible for the production 
of about half of the observed abundances of elements heav- 
ier than iron. Although many sites have been proposed for 
the r-process, the neutrino energized wind above the proto- 
neutron star in core-collapse supernovae remains a favorite 
dWooslev et al.| [T994). Whatever the environment, however, 
it can be shown that the solar-system r-process abundances 
are well reproduced by beta-decay flow in a system that is in 
approximate (n,7) ^ (7,n) equilibrium. Hence, the relative 
abundances of isotopes of a given element are determined by 
the revised nuclear Saha equation ( l24l i 



n{A) 
n(A+l) 



^ >nA ; Ga+i 



(57) 



(58) 



where Q„ is the neutron-capture Q-value for isotope '*Z (or 
equivalently the neutron separation energy for the nucleus 
^^^Z) and n{A) represents the abundance of an isotope ^Z. 
This equation defines a sharp peak in abundances for one (or 
a few) isotopes within an isotopic chain. The flow of beta de- 
cays along these peak isotopes is then known as the r-process 
path. 

T he location of the r-pr ocess path peak is roughly identi- 
fied (iBurbidge et al.|[T957h by the condition that neutron cap- 
ture ceases to be efficient once n{A + \ )/n{A) ^ 1. Taking the 
logarithm of Eq. (ISTl l and inserting the numerical terms, the 
r-process path can be identified by the following relation: 

OA 

J path 

/ 3 kT n 
2.30 35 .68 + -log ( )-log(^) + log(l +/?) 

\ 2 MeV cm ^ 

For a typical r-process temperature of ^ 1, the require- 
ment that the r-process path reproduce the observed abun- 
dance peaks at A = 80, 130, and 195 implies that the r-process 
path waiting point in beta flow halt near the neutron closed- 
shell nuclei '^"Zn, '^"^Cd and ''^^Tm. For a neutron density 
sufficiently high (n„ ^ 10^" cm"^) so that the neutron-capture 
rates exceed the beta-decay rates for these isotopes the peak 
abundances along the r-process path must be for isotopes 
with Q„ ^3 MeV, and thus (Q„/kT)patb ^ 30. For this ra- 
tio, the correction factors in Eq. 02] | are negligible. This con- 
straint on Q„, however, concerns the conditions near "freeze- 
out" when the final neutrons are exhausted at the end of the 
r-process. At this point, the system falls out of NSE and the 
r-process path decays back to the line of stable isotopes. 

From Eq. ( |58] ) we deduce that Q/kT ^ 1 along the path re- 
quires a neutron density of ^ 10^"' cm"^. Although such con- 
ditions do not occur at freezeout, they may occur earlier in 
the r-proce ss. For example , in the neutrino driven wind mod- 
els of Wooslev et al.l (11994) . the r-process conditions begin 
with the neutron density at n„ sa 10^^ cm~^ and a temperature 
of Tg ~ 2. The density is also much higher (> 10^^ cm"^) 
when the material is first ejected from the proto-neutron star. 



Such conditions may also be achiev ed for an r-process which 
occurs during neutron- star mergers dFreiburghaus et al.ll 19991: 
iRosswog et al.|[T999ll2000l) . 

As the correction factors become larger, the effect on the 
r-process may be to slightly increase the time it takes for the 
r-process to build up substantial abundances of the heaviest 
nuclei. This is because the faster photodisintegration rates 
will drive the r-process path slightly closer to the line of beta 
stability. 

As an illustration. Figure |7] shows the reverse ratio cor- 
rection factors for some neutron-capture reactions relevant 
to the peaks in the r-process abundance distribution near 
A = 80 and A = 130. These plots we re generated from 
Eq. (EB using the REACLIB compilation dCvburt et al .120101) 
for [Na{<Jv)(T<))] and reaction Q-valnes. Correction factors 
for the ^°Zn(n,7)**'Zn (Q = 2.151 MeV) reaction and the 
'^°Cd(n,7)''"Cd (Q = 2.028 MeV) reaction are relevant to 
the r-process path near freezeout and are quite negligible 
at the termination of the r-process for T() ^ 1-2. Correc- 
tion factors for the **'^Zn(n,7)'*^Zn (g = 0.240 MeV) reaction 
and the '^''Cd(n,7)'"Cd (Q = -1.5 MeV) reaction are rele- 
vant to the r-process early on when the neutron density can 
be very high. In this case, the correction factors are ~ 5% 
for 79 ~ 1 - 2. For comparison, the dashed line on the plot 
for the '■"^Cd(«,7)'^'Cd reverse reaction is from the leading 
term in the analytic nonresonant neutron-capture expansion 
of Eq. (|52] |. This shows that the simple analytic expression 
accounts for most of the correction factor. Note also, that the 
'""'Cd(n,7)'-'^Cd involves a negative Q-value. This causes the 
correction factor to be ~ 5% even down to low temperatures 
Tg < 0.5. 

3.2. rp-Process 

The rp-process (IWallace & Woosleyl 119811 : ISchatz. et alJ 
Il998h refers to a sequence of rapid proton and alpha cap- 
tures in explosive hot hydrogen-burning environments. The 
most energetic such environments occur on accreting neutron 
stars and are believed to be the source of observed type-I X- 
ray bursts. Energetic events are also thought to occur during 
accretion onto the magnetic polar caps of X-ray pulsars. In 
such environments rapid proton captures and beta decays oc- 
cur among proton-rich isotopes until further proton capture is 
inhibited because such capture would lead to an unbound nu- 
cleus. At this point, the process must wait until the last bound 
nucleus can beta decay so that further proton captures can oc- 
cur. This pro cess continues unti l it is terminated for nuclei 
with A - 100 ("Schat z et al.ll200T ). 

Proton captures in the r/7-process are sufficiently rapid that 
(/?,7) ('j,p) equilibrium can be assumed. Hence, the rela- 
tive abundances of isotopes along a given proton-capture path 
are again determined by the revised nuclear Saha equation 
(EH 



n(Z,A) 



n(Z+l,A+l) 



■= l+R] 



1 / l-K^LklX^^^ GzaG„ 



«A ) 



-Qp/kT 



(59) 

where now Qp is the proton-capture 2" value for isotope ^Z 
(or equivalently the proton separation energy for the nucleus 
^^'Z+ 1) and n(Z,A) represents the abundance of an isotope 
^Z. Just as in the r-process, this equation defines a peak in 
abundances along an isotonic chain. In this process, proton- 
capture rates along the r/9-process tend to be dominated by 
resonances in the captured nuclei. As noted above, the quan- 
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turn corrections will be largest for nuclei with low proton 
binding energies and also low lying resonances. 

Recently, the reactions which most significantly affect the 
r/?-process models have been reviewed (Parikhetal. 200^. 
Correction factors for some of the most important proton- 
capture reactions identified in Parikh et al. (2009) are illus- 
trated in Figure [8] These plots were generated from Eg . OTT l 
using the REACLIB compilation dCvburt et alj |20T0 ) for 
[Na{(Jv){Tc))] and reaction 2-values. Shown are correction 
factors for the ^'^S{p,j)^^Cl (Q = 0.296 MeV), 42Ti(/?,7)'*3v 
(Q = 0.087 MeV), '^Crip^-ff'^Mn {Q = 0.467 MeV) and 
^Ge{p,-f)^^As (Q = 0.169 MeV) reactions. Note that the 
2-values listed here are fr om the REACLIB t able and dif- 
fer from those estimated in iPariMi et al.l (l2009h . For typical 
rp-pmcess conditions (Tg ^ 1 -2.5) the correction factors are 
^0-3%. The effect of including these corrections will be to 
increase the photodissociation rates and therefore to shift the 
path slightly closer to stability. This could slightly modify the 
time-scale and energy release in X-ray burst models. 

3.3. Core and Explosive Neon, Oxygen, and Silicon Burning 

Once the core (or shell-burning) temperature in a massive 
star exceeds a temperature ^1, photonuclear reactions be- 
gin to dissociate nuclei into free protons, neutrons, and alpha 
particles plus heavy nuclei. This leads to photonuclear rear- 
rangement as the free nucleons and alpha particles are cap- 
tured on the remaining nuclei, eventually leading to the build 
up of an iron core. This process begins with core oxygen 
burning (Tg ^ 1.5-2.7) and culminates with core or explo- 
sive quasi-equilibrium silicon burning which occurs at tem- 
peratures of up to ^ 5. 

This process can be described by the NSE equation. For 
the most part, however, this process involves photonuclear 
rearrangement of nuclei near the line of beta stability with 
rather large nucleon and alpha separation energies, 2^5-10 
MeV. Hence, even up to ^ 5, the major products of these 
advanced burning stages involve Q/kT > 10 and the correc- 
tion factors can be neglected. Nevertheless, there are a few 
minor reactions during the photonuclear rearrangement for 
which Q ^ 1-2 MeV. Some examples of reactions which 
are slightly affected during core and explosive silicon burn- 
ing are shown in Figure |9] These include the reverse rates for 
the ^"^Mgip , 7)25 Al (g = 2 .27 1 MeV), ^^Si(p, -ff^P (g = 2 .747 
MeV), and ^^Ar(p,7)"Cl (Q = 1.857 MeV) reactions. Even 
for these nuclei, however, the maximum correction is < 1%. 
Hence, for the most part these corrections can be neglected 
during core and explosive thermonuclear burning in massive 
stars. 

3.4. ^-Process 

Perhaps, what may seem as the most obvious applica- 
tion of the d educed corrections w ould be to the 7-pr ocess 
(IWoosley & Howardl978;.Hayakawaet al.ll2004l200Q) . The 
nucleosynthetic origin of the isotopes th at lie on the proton- 
rich side of stabihty lik ely requires (IWooslev & HowardI 
1 1 978t iHavakawa et al.l [20b4, 2006) the onset of photonuclear 
reactions with Tg ~ 2-3. In this scenario, the nucleosynthe- 
sis of proton-rich nuclei is achieved by photo-disintegration 
reactions on pre-existing s- or r-process nuclei. The nucle- 
osynthesis of all relevant nuclei in this process, however, in- 
volves a path on the proton-rich side of stability for which the 
proton binding energies are of order Qp ^2.6 MeV, imply- 
ing Q/kT > 5. Hence, the correction factors deduced here 



have little effect on the production of proton-rich nuclei in the 
7-process and can be neglected. 



3.5. Big Bang Nucleosynthesis 

For the most part, standard big bang nucleosynthesis is 
also unaffected by the corrections derived here. To exam- 
ine this we have calculated big bang abundances with and 
without the correction for quantum statistics. The temper- 
ature evolution of the light element abundances is summa- 
rized in Fig ure [TOl To gen erate this plot, we used the net- 
work c ode in lKawanol d 1 992h and Smith, Kawano, & Malaiie^ 
( Il993h with reacti on rates from Descouvemont et a}. (20Q# 
and ICyburt & Da vids (20 081). The adopted neutron life time 
is r„ = 881.9 s (iMathewsTKajino. & Shimal l2005h and the 
baryon to photon ratio is fixed to 77 = 6.2 x 10~^° based upon 
the WMAP seven y ear data for the ACDM-nSZ-nlens model 
dLarson et al.ll2010l) . 

Because the light-element reactions involve relatively high 
2-values and most of the nucleosynthesis involves relatively 
low temperatures Tg < 1, there is almost no effect on the fi- 
nal calculated abundances. Even at the epoch during which 
the weak reactions decouple (Tg ^ 10), the slightly modified 
NSE abundances (see insert on graph) would have a negligi- 
ble effect on the details of neutrino last scattering. Hence, 
one is justified in ignoring these corrections for standard big 
bang nucleosynthesis. It is possible, however, that some sce- 
narios of inhomogeneous big bang nucleosynthesis could in- 
volve proton-rich or neutron-rich nuclei with lo w Q-values 
dMalanev & Fowler |[T989t iKajino & Bovd |[T990h . 

4. CONCLUSIONS 

Since photons are massless, they are sensitive to the ef- 
fects of quantum statistics. As a result, at high tempera- 
ture some of the reverse reaction rates with low Q-valnes 
can be modified from the tabulated values based upon the 
approximation of classical Maxwell-Boltzmann statistics for 
the photons. Moreover, in any environment the quantum ef- 
fects always speed up the photodisintegration rate because a 
Planck distribution places more photons at low energies than a 
Maxwell-Boltzmann distribution of the same temperature. As 
a result, those nuclei with loosely bound nucleons lose them 
faster than predicted by the classical Maxwell-Boltzmann dis- 
tribution. In this paper, we have derived analytic correction 
terms for the effects of the quantum statistical distribution of 
photons on tabulated thermonuclear photodisintegration rates. 
Usually, this modification is small because (Q/kT) ^ 1 so 
that there is little difference between a Planck and a Maxwell- 
Boltzmann distribution. The effect is largest for environments 
for which synthesized nuclei have Q/kT ^ 1, however, the 
correction is still small compared to the uncertainty in the es- 
timated thermonuclear reaction rates for such nuclei. 

We have analyzed possible effects of these corrections in a 
variety of astrophysical environments including the neutron- 
capture r-process, the hot hydrogen burning rp-process, core 
or explosive silicon burning, the photonuclear 7-process and 
big bang nucleosynthesis. In general these corrections have 
little effect except, perhaps in the case of the rp-process for 
those reactions near the proton drip line waiting points, or in 
the early stages of the r-process when the neutron density is 
high enough to drive the r-process path to nuclei with low 
2-values even at high temperature. 
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Figure 2. First few con'ection coe fficie nts for cliarged-pailicle reactions as functions of Eg / kT for tli e ca se of a constant S-factor. Tlie dashed lines show the 
approximate values based upon Eq. (43) and the solid lines show an exact numerical integration of Eq. ^28). 
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Figure 3. Reverse rate correct ion factor R for charged-particle reactions as a function of Q/kT for various values of Eq/Q. Tiie soli d line s in this figure are from 
a numerical integration of Eq. I14t for the case of a constant S factor. The dashed lines show the result of the series expansion in Eq. (46) with only the first three 
terms included. 




Figure 4. Reverse reaction rate coirection factor R for the hydrogen-b urni ng charged-particle reactions ^Be(p,7)'*B (Q = 0. 135 MeV) and ''He(a, 7)'Be (Q = 
1.587 MeV) as a function of Ti). These plots were generated from Eq. )3U using the REACLIB compilation I Cybur t et a l. 2010) for [Nj^{av)(Tg)] and reaction 
Q-values. On each plot, the solid line shows the application of Eq. J3U . the dashed line shows the application of Eqs. )44K while the dotted line shows the results 
obtained from keeping only the first term in the series in Eq. j46). For the ^Be(p, 7)'*B a constant 5-factor was assumed, while the dashed line for the^He(Q:. 7)'Be 
reverse reaction is based upon (S(0)/5(0),5'(0)/25(0)) = (-0.953 MeV"', 0.835 MeV-^) from lCaughlan & Fowleil (198^) . 
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Figure 5. Reverse reaction rate correction factor ff as a function of Q/kT for nonresonant neutron-capture reactions with the assum ption of a slowly varying S 
factor The solid line is from an exact numerical integration of Eq. )14t . The dashed line is from the analytic expression given in Eq. )52t truncated after the first 
three terms. 
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Figure 6. Reverse ratio correction fac tor R for the resonant hydrogen-burning charged-particle reaction '-C(p.7)''N (Q = 1.944 MeV) as a function of Tg. The 
solid hne was generated from Eq. BU using the REACLIB compilation iCy burt et al. 2010) for [Na {crv)(Tg)] and reaction 2-values. The dashed line shows the 
result from application of the simple single resonance correction in Eq. i56\ for a resonance energy of = 0.434 MeV. 
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Figure 7. Correction factors R for the reverse rat e of some neutron-capture reactions r elevant to the peaks in the r-process abundance distribution near A = 80 
and A = 130. These plots were generated from Eq. 1311 using the REACLIB compilation jCvburt et alJ2O10l) for [NA{cn')(kT)] and reaction Q-values. Correction 
factors for the *''Zn(H,7)*'Zn {Q = 2.151 MeV) and '^''Cd(«,7)'^^'Cd {Q = 2.028 MeV) reverse reactions are relevant to the r-process path near freezeout and are 
quite negligible at the termination of the r-process for ~ 1 - 2. Correction factors for the ^^Zn(n. 7)^'Zn (Q = 0.240 MeV) reaction and the '^*Cd(«, 7)'^'Cd 
{Q = -1.5 MeV) reaction are relevant to the r-process early on when the neutron density and temperature can be very high. In this case, the correction factors are 
up to ~ 5% for ~ 1—2. For comparison, the dashed line on the plot for the ^^^Cd(n. 7)'^' Cd reverse reaction is from the leading term in the analytic expansion 
ofEq. (52j. 
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Figure 8. Correction factors R for the rev erse rates of some proton-capture r eactions important f or the rp-process in models forZ-ray bursts IParikh et alJ2009h . 
These plots were generated from Eq. 13 It using the REACLIB compilation jCvburt et alJI2O10l) for [W^ (crv) (Tg )] and reaction Q- values. For typical rp-process 
conditions {Tg ~ 1-2.5), the correction factors are ~ 0-3%. 
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Figure 9. Correction factors R for the revers e rat es of some proton-capture reactions which occur durin g core or explosive oxygen or silicon burning in massive 
stars. These plots were generated from Eq. j3U using the REACLIB compilation jCvburt et al.ll2O10l) for [NA(uv)(Tq)] and reaction Q-values. Even for the 
maximum silicon burning conditions (Tg ~ 5), the coirection factors are < 1%. 




Figure 10. Calculated BBN abundances as a function of cosmic temperature with (solid lines) and without (dotted lines) the quantum coiTections. The insert 
shows an expanded view of the deuterium abundance in the interval from = 100 to 10. 



